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ABSTRACT 



In the heterogeneous flow of finely divided particles and gas 
in a rocket nozzle, the assumptions that the flow enters the nozzle 
with no lags, that the radiative transport of energy between particles 
is negligible, and that there are no temperature gradients within the 
particles, are commonly made. In this paper, each of these as- 
sumptions is investigated in detail, and they are shown to be reason- 
able for typical rocket nozzles. 



I. INTRODUCTION 



In recent years, it has become common to add metals to the 
propellant in rocket engines to increase the chamber temperature. 
Unfortunately, most of these additives form some sort of a solid 
particle in the exhaust, and consequently can cause a depreciation in 
the thrust and specific impulse from the values based on gaseous 
products. That is, instead of an expected 20 per cent increase in im- 
pulse, say, an increase of 12 per cent is noted. 

A very common mixture is a propellant of ammonium per- 
chlorate oxidizer and a poly sulfide fuel with aluminum added. For 
the numerical work in this thesis, this propellant is assumed with 
pure aluminum added to it. 

The amount of metal that can be added is limited, since there 
is a requirement for the fuel to act as a binder for the oxidizer and 
metal. Usually, the designer will add as much metal as he can; thus 
a method is needed whereby he can correlate test engine performance 
with the predicted performance. 

Such a method is presented by Rannie in reference 1. In ref- 
erence 1, Rannie used a one -dimensional hydraulic approximation, 
which is known to be valid for homogeneous flow, for the heterogeneous 
flow of gas and solid particles. In addition, he assumed (i) that the 
flow entered the convergent section of the nozzle with no lag between 
the ga3 and solid particles; (ii) that the radiative transport of energy 
between particles was negligible; and (ill) that the thermal conductivity 
of the solid particles was sufficiently high compared to that of the gas 
that there were no temperature gradients within the particles. 
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As sumption (i) is Investigated in Part II by solving the equa- 
tions of motion in a cylindrical combustion chamber and determining 
the velocity lags at the entrance to the nozzle. Axially symmetric 
flow is assumed in the chamber, and the mixture is approximated by 
a heterogeneous mixture of gas and solid spherical particles. A 
constant source distribution of gas and particles is asstlmed along 
the cylindrical wall, and viscous effects in the main flow are ne- 
glected. Turbulence and temperature gradients in the chamber are 
also neglected. It is shown that the particle lag at the entrance to the 
nozzle is indeed very small compared to the gas velocity, and thus 
the assumption that the lag is zero is a valid one. 

The magnitude of the effect of radiative transport of energy 
between particles is estimated by adding a radiation term to the equa- 
tion for the heat balance of a particle and investigating the size of 
this additional term. The ratio of the absorption cross section to the 
geometric cross section is a necessary parameter in this term, and 
a value is not available for aluminum oxide particles of these small 
sizes and high temperatures. However, for any reasonable value of 
this parameter, the radiative transport of energy is shown to be neg- 
ligible compared to the heat transfer by convection and conduction. 

The temperature distribution within a spherical particle is 
found for unsteady heat flux over the surface of the particle in terms 
of the Laplace transform. By expansion of the transform, the solu- 
tion for uniform temperature within the particle is recovered and the 
magnitude of the correction term is determined. The correction term 
is shown to be negligible for practical cases. 
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II. WO- DIMENSIONAL ANALYSIS OF THE FLOW 
IN A COMBUSTION CHAMBER 



Following the notation of reference 1, let oC be the mas 3 
fraction of particles in the heterogeneous mixture, fy the gas den- 
sity, and fe the density of the solid particles. Then the density of 

the mixture is given by 

I I — <*■ 

r * f } (5 (1) 

Let A be the cross-sectional area of the cylindrical combus- 
tion chamber and m the constant mass flow rate of the mixture. Ap- 
plying the condition that the net flow rate into any volume must be 
zero, the continuity equationofor the solids and gases separately are: 



fr [0-*) N 




( 2 ) 



_cL 

3a 





( 3 ) 



where U.^ and are the axial velocity components of the gas and 

solids respectively, and Vo and Vj are the radial velocity com- 
ponents of the gas and solids. 

By considering the momentum transport through a small ele- 
ment and neglecting friction forces, the momentum equation in the 
axial direction can be written 



JL 

3 x 









+ iLf * pilyV s K] = -^£-r 
L -I 3 a. 



( 4 ) 



4- 



and the radial momentum equation is 

’ rp [('-’ / fl/ ? l r ] -"4^ + ? 

Approximating the particles by spheres of radius cl , and 
assuming that a modified Stokes drag law holds with a correction fac- 
tor fj to allow for deviations from Stokes' flow, the force balance 
in the axial direction on a particle is 

f *^>37-+^]= fc-rr ^(U r U,)j- - (6) 



The steady 6tate drag law is assumed valid. JLc is the viscosity of 
the gas and is constant since the temperature is assumed constant, 
and the force due to the pressure gradient is added. The coefficient 
fj is a function of the Reynolds number and the Mach number of the 
gas relative to the particle, and is equal to 1. 0 in the Stokes flow re- 
gime. 

A similar equation for the radial force on a particle is 



f iro. 3 P s [u,^i + v; <^j = -i- ir a 3 ^ < 7) 

Equations 6 and 7 are expressed more conveniently in the fol- 
lowing form: 









^ J^U-s ~ ~ ~ "t" ^4 



/* 






+ 



?S 




(8a) 
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and 



fl ( UfAli. + Vr A^ 4- -L- i^p2-) 

^ Z 4 ( 3 x dr f s dr ( 



( 8 b) 



Letting U s —U^ f (U.j - \JL^j and \/ 5 — in equation 

3 and subtracting it from equation 2 









( 9 ) 



Introducing the following non-dimensional variables <£" 77 - • 
{ f= -dd- , ^ = _t_ , where V 0 is the radial ve- 

locity component of the gas at | just outside of the boundary lay- 
er and b is the radius of the burning surface. To simplify the 
equations, also let 

e = 2. ^ Mo fj 

^ yU b 

Introducing these in the above equations, we obtain for the combined 
continuity equation 




( 10 ) 



and for the momentum equations 



if[( 1 **] 7 + ^[(' V] + 41" f 7 

+41 <f +*«?]- --&§*• 



( 11 ) 



(12) 



■ 6 - 



'Mrr^] = 



'U 



v* V. 



and for the force balance equations, 



V +s " 



+. vy ^ , J_ yp' 

** J M ?x ^ j 



(13) 



V’'“ e j>T^ 



<6-0- ♦-» 



i£. 



^ fX^7 



(14) 



With the help of equations 2 and 3, equations 1 1 and 1 2 can be 
further simplified by rearranging so the dependent variables are , 
^ » 4$- 4^ . and f s - fa . Thus equations 11 and 12 become 

♦ V-^l + ( ^' 

F< V^s" t + *3 

+ r^\ - zX 



First Approximation, Zero Velocity Lags. 

When the velocity of the particles approaches that of the gas, 
(i. e. , as CL ° ), 6 — > o and <p t — >. and 4^. 
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Then 

Cp s = = r and Ts= tg = ^ 

where the subscript zero indicates no lag. 

Then equation 10 reduces to 

4. Aik . Jk_ _ o 
2 ) f r £>7 V, “ 

and equations 15 and 16 reduce to 




(17) 



(18) 



(19) 




( 20 ) 



It is a simple matte r to verify that the solution to these equa- 
tions, satisfying the boundary conditions =0 at f =■ o 
and H'o— -/ at ''l ~ 1 , is 



4>o- Zf 
% = -7 



_f_ = -zr -i?' 



+ Cc n £X • 



(21) 



Second Approximation, Small Velocity Lags. 

We shall assume for the solution that the variable may be 
expanded in powers of 6 ; thus 
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= + e fy, 






‘f’s, 


+ 




- v » + 4 r J, 


+ 


n 


- + 6 V ( , 


+ 




= + 6 °C, 


■t 


r 


= r e +« f, 


+■ 




J 



The criterion for these expansions to be valid is that 



(22) 



<£o fo 




(23) 



and similar inequalities for the other variables. Here 6 is merely 
an indicator of a small quantity. Substituting these expansions into 
equations 13 and 14, and equating coefficients of £ , 






- <*>. 



2>f, 






+ 2> 7 






and 









+ ^ a 



2 Vo 

1 



- f« 






The density of the misrture, f’e, , is very much less than the density 
(° s of the solid material, so those terms can be neglected that con- 
tain » therefore 











(24) 



and similarly 
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1% 









(25) 



Now defining the first two terms of equation 10 as p , then 






~1 



(26) 



Equation 2 can be written as 

F ^ L* p V i 

Then differentiating this by parts, 

r a 



=- o 



F - °C 



■iH 



= o 



or 



— o 



^'-<9 -f’+jif -^isf 

Also comparing equations 26 and 10 we see that 



(27) 



F - 









(28) 



From the first approximation, we may evaluate equations 24 
and 25 as 

tj,- %, = >? 



( 29 ) 
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Substituting these into equation 28, the solution of F is 

F = 



( 30 ) 



where use was made of equation 22 and only carried to order epsilon. 

Putting equation 30 back into equation 27 and incorporating 
equation 22, and rearranging. 



- ^oLi 







(31) 



Particular solutions for equation 31 are 

, = 3 °< 0 ( I — c(F) f -h A 



(32a) 



eC, = -L <*.*(/ -6* J f 



8 



(32b) 



The boundary condition for is chosen so the mass flow 

rates of gas and solids are the same as with no lag. . Since 



(33) 



As. = e 0 & & 

A L 



v o 

the conditions that ^ ^//\ and rn \ are the same as with no 
lag, are 



% 

SLi+f^^rf = 
■«t. f. ' 



(34) 



From equation 1, upon substitution of equation 22, and .dropping 

Z- 



terms of order ^ 



9 
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K ( 

I 



( 35 ) 



Substituting equations 35 and 29 into equation 34, we find that 

-oCo(l-^-o) ->1=1 ( 36 ) 

satisfies the condition that the last term of equation 34 vanish. 

Clearly the solution 32a is not applicable, and applying condition 36 
to equation 32b, 

oC i - ~ L oLo ( y f -L] (37) 

The boundary layer and friction forces along the wall are ne- 
glected and the mass fraction of particles to the mixture is assumed 
to be constant along the length of the chamber since the propellant is 
presumably uniform in mixture. 

This solution does become infinite at the centerline, so the 
solution is valid for the whole chamber except on a line down the cen- 
ter of the chamber. This singularity does not cause any difficulty 
because the mas 3 flow rate in the axial direction is finite, as will be 
shown below. 

Now substitute equation 2 2 into equation 15, and put -P^fo* > 
then equate the coefficients of & , dropping known zero derivatives 

and ); the result is 

* ^ ^ % 4r^- 

dy 




Jr. j p / n f y /f \ c- ^ ° — — L_ fbfr- 

oio?e[^ l 4 ^,) MS" ^ 
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and upon substituting equations 21, 29, and 35, 









? + ^ ^ - «• - ij& 



H 

b Y 









Combining the second and third terms and putting in the value (equa- 
tion 37) of °(~i , we obtain 






Similarly, equation 16 becomes 



7 






(38) 



+x -^ 7 !- r 77-7^ (39) 

Let us assume a solution of equation 39 to be 

^fi = ~ 3 1+1^ (40) 

Substituting equation 35 into the continuity equation for the gas (equa- 
tion 2) leads to 

^9_l _ _ -L )(l3A 

(41) 



Mj. _ _ -L jjl % j} 



Upon substitution of into equation 41 and integration, 

fy, — (> i~3°Cc ? (42) 

where the constant of integration disappears since 4^, must be 
zero at f — ° 



From equation 29 we find 






and 



(43) 
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To verify the solutions, substitute the values for , 4* g./ > 

tg-j » and 4 S , (equations 40, 42, and 43) into the continuity equa- 
tions (equations 2 and 3). The mass flow rate of gas for a cylindrical 
chamber of radius b , with no lags, is 

— fr ^(i fo Vc 4> e (44) 

and v/ith velocity lags, is 

= X-fl r £ (l Fo Vc to j ( I +* e ^l X )] ^ 7 (45) 

J 6 



Substituting from equation 42 and integrating, equation 45 becomes 
m u = (46) 

The mass flow rate of the particles without a velocity lag is 

= fr b oLo^o^o <P 0 (47) 

and with a velocity lag. 




o_rr b c^PoVo 





(48) 



Evaluating the integral 

^ (49) 

showing that the mass flow rate of the particles does not change when 
velocity lags are introduced. Hence, the assumed solution and con- 
stants are consistent with mass conservation. 

To obtain some numerical results, a specific chamber must 
be chosen. The following values are consistent with the numerical 
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example given in reference 1. 

b =o, 3 io (c ft 

JUL = / ,5"* lc° 

f l - Zlo ly /f c * 

T = booo °P 



I 

CL = Z K 1 o c 



r - o.s~ '"/sec. 

t? 



~f a ■=. I ooo ^/n 



< = V 



m . 



The chamber length is chosen a3 twice the diameter, and Stokes 
flow will be assumed f l.o) . Then 

Vo = -| 2n tysec. 4 , o= 7 lS>%ec 

6 = 3.H3 * 'o' 3 % - -/ 9.7 ^jec, 



The magnitude of the velocity lags in relation to the no-lag 
flow velocity is obtained from equations 29 and 21 

(50) 

which verifies that the lags are indeed small. 

The uniformity of the flow across the exit plane was assumed 
in the analysis in reference 1, and Table 1 is presented as a verifica- 
tion of this assumption. Values of <= ^fy1p 0 , e <f o , and 

were calculated as functions of y , the non-dimensional 

radius. 

The Reynolds number of the gas relative to the particles was 
computed at the exit plane using the particle diameter as the charac- 
teristic length and found to be 6. 5 X 10 The Mach number was 

-4 

also found to be 1. 36 x 10 . Combining these two values, a value 

_3 

for M/Re of 2. 09 x 10 is found, which indicates that the particles 
are in the Stokes flow regime, and thus the correction factor fy is 
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TABEE 1. 



\ 








4 


<L 0 


. 1 


-7. 35 x 10* 4 


-3.69 x 10* 4 


2. 63 X 10* 2 


. 2 


-6. 10 x 10* 4 


-2. 43 x 10* 4 


1. 73 X 10*‘ 


. 3 


-4. 75 x 10* 4 


-1. 55 x 10* 4 


1. 28 X 10* 


. 4 


-4. 45 x 10* 4 


-9. 15 X 10* 5 


9. 25 x 10* 3 


. 5 


-4. 07 x 10* 4 


-4. 21 x 10* 5 


6. 49 X 10* 3 


.6 


-3. 69 x 10* 4 


r 

-2. 29 x 10*° 


4. 24 x 10* 3 


. 7 


-3. 35 x 10* 4 


+3. 14 x 10* 5 


2. 35 x 10* 3 


. 8 


-3. 05 x 10* 4 


6. 15 x 10*‘ 


6. 95 x 10* 4 


.9 


- 2. 8 1 X 10 


8. 41 x 10* L 


-7. 58 x 10*' 


1. 0 


-2. 55 x 10* 4 


1. 10 x 10* 4 


-2. 06 x 10* 



approximately equal to one, as taken above. Reference 1 gives a dis- 
cussion on the calculation of the parameter rj . Using equation 66 
and equation 67 of reference 1, 

C =. I ,e> •+• H-SS = l.ooU (51) 

Thus, from the above analysis, the assumption made in refer- 
ence 1, i. e. , that at the entrance to the nozzle the flow was uniform 
and without lags, is a reasonable assumption. 
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UI. RADIATION BETWEEN PARTICLES 



For a one-dimen3ionaI steady flow of a homogeneous fluid 
without heat addition, the energy equation, in terms of internal energy 
per unit mass, is 

f (\:V- ^ f ^ = Centre, (52) 

and for the heterogeneous mixture of gas and particles is 



(^t^y*e4Ll s (e 5 r r u s 

3 



C on ST . 



(53) 



where oC is the mass fraction of particles (subscript S ). This 
equation is general and independent of force interactions and heat 
transfer processes between the particles and the gas. 

If radiative transfer of heat between particles is allowed, an 
additional term must be added on the left hand side of the equation, 
since the equation does not take care of interactions and heat transfer 
processes between particles. The energy transfer between particles 
due to radiation is 



where 




(54) 



(55) 



may be regarded as an effective conductivity due to radiant energy 
transfer. The coefficient CT^ is the Stefan- Boltzmann constant, 
and Xu is the Roeseland mean free path. Then — f\ 

must be added to the left hand side of equation 53 to account for radi- 



ative transfer. 



Equations 54 and 55 are applicable only when -Zjz is very 



small compared to the length scale appropriate to the problem, for 
instance the nozzle throat radius, at all important wavelengths of 
radiation. 

The variable mass fraction c-(_ in equation 53 can be elimi- 
nated with 



obtained from equations 2 and 3. In cases of interest, the volume oc- 
cupied by the particles is negligible when compared to the volume of 
the gas, and the assumption 



prosdmation is valid over a wide range, but should be checked for any 
particular problem. 

Substituting equations 56 and 57 into equation 53, the energy 
equation becomes 




(56) 




(57) 



3 

is valid for a particle material density of 210 lb/ft as we have chosen. 

3 

Since the gas density will seldom be greater than 0. 4 lb/ft , this ap- 





( 58 ) 



where ^ c is the enthalpy In the chamber. 



For a perfect gas with constant specific heats 
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and we now introduce a specific heat, , for the mixture with a 

particle mass fraction o^ 0 , 

Cp^ =• C l—cCo)£p + ct-o £ (60) 

and modify the gas constant to 

R 0 = (i -**) R (6i) 

The right hand side of equation 58 is where T c is 

the stagnation temperature of the mixture. Then equation 58 may be 
rewritten as 



(l-UdpTi-ttoCT, t (i-^)dJs = tpcJ d (62) 

11 ^ ^ 1 ' 4 * 1 

where equations 56 and 59 were used. 

Now writing T s as Tcj + l s — and Lii as 
and combining the terms of equation 62, 



(f +*.4T$-T ? }+ = Cp„T 0 



(63) 



Neglecting radiation, the heat balance for a sphere of radius 
and temperature is 

^-fra. f^C U s ( Ts -T 9 ) ‘i < 64 > 

where is the conductivity of the gas surroimding the sphere and 

is a correction factor that takes into account the influence of con- 
vection. For the Stokes regime of flow, =. ),£> , and it decreases 



at higher Reynolds numbers. 
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Rearranging equation 64, 

T % -T s = 

t ■ k t 3T 



( 65 ) 



To express equation 65 in dimensionless form, the following 
notation from reference 1 is introduced 



U- = \TrTtT t 

r- Tc r 

X - r r f 

. _ x cL x r<\ r ^fh 
" 1 r r /fe 

giving 

T,-T S 

where 




( 66 ) 



(67) 



( 68 ) 



has a numerical value not much different from unity. If the particle 
radius is extremely small, £ approaches zero and since the other 
factors on the right hand side of equation 67 are of order unity, 

and the temperature lag is zero. The velocity lag is 
also proportional to e so it becomes zero as well. As increases, 
e becomes larger, and the principle of the perturbation procedure in 
reference 1 was to expand ^ , I s ' ^ in power 



series in £ . 



CT( ro 
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V/ith radiative transfer, there is an additional term in the lie at 
balance for the particle. Suppose a particular particle is located at 
X = \ 0 • The radiative transport of energy in the X - direction per 

unit area normal to X in the space Surrounding the particle and close 
to the particle is 




Hence the net radiative transport of energy to the sphere is 



IV 




and the heat balance for the particle becomes 




which replaces equation 64. Rearranging and introducing the dimen- 
sionless variables from equation 66 









(70) 



The radiative mean free path ^ £ is given by 






n 'ft'- cCr 



( 71 ) 
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whcre n is the number of particles per unit volume and GL is the 
ratio cf absorption cross section to geometrical cross section. From 
the identity 



n±n-«?fs = 



and the equation of state for the gas 

-± 



n 



= 3 7 % 



c< 



Y / -o- \ZT« f s a, Y ^ /-o( 









Substituting into equation 7 1 above 



= 4- ^ 6*- — 









r 



3 fa. <2i - 

and hence from equation 55 becomes 

-l ^ r %-/ 

^ ^ S Gc ^ 

Substituting for ^ into equation 70 



(72) 



(73) 



(74) 



(75) 



Y- d . 32. p GlSjj 

- r *lrt nHs ji + 3 1 



c a. _l_ d 






!76 > 



where Q is assumed to have small variation with and can be 
t alcen out of the differential. Since fy/g, ~ v<-i , equation 76 

can be written as 



V** if 



(77) 



f = 



<r Ifi? a ( 



r fe WoT c rv dt 



where 
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Sincc the factors multiplying j- in equation 77 can be con- 
sidered to be of order unity, as is the first term in the square brackets, 
a value of of order unity implies that radiative transpoi't is as ef- 

fective as the heat transfer by convection and conduction. In the exam* 

~ ~ q / 

pie of reference 1, 3*1°/ & ; hence, -fy is extremely small 

even if Q_ is quite small. Hence the influence of radiative transport 
on particle lags is negligible. 

The zero lag solution is also affected by radiative transport of 
energy. The energy equation 63 becomes 




kir ( dJjz . 



Upon substitution of equation £5, we find that the condition for radiation 
to have negligible effect on the zero lag solution is 

i (78) 

te/fca ^ 

since the non-dimensional terms are all approximately of order one. 

The concept of an effective conductivity has meaning only if 
r r ^ / . Tor chemical rockets ^ R a r<- ) = 1° 

Hence the radiative term gives negligible contribution. The ratio 
is 

^ 3 ^ d RT1 J- 

In the e;cample ■= 2 s I o S A / 

Zmi — US' • Even if is small, / . 

c\ / cr 



, acid 
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IV. EFFECT OF PARTICLE CONDUCTIVITY 

The effect of finite particle conductivity is investigated by 
solving the heat conduction equation for the particle with an appropri- 
ate boundary condition, then expanding that solution to investigate the 
size of the correction term. 

The heat conduction equation for a spherical, uniform, solid 
particle with spherically symmetric temperature distribution is 



IX = r W= crJ^f r ( r^) 

where 0“ is the diffusivity of the particle material. 
The initial condition is 

T = at t = o 

and the boundary condition at the particle surface is 




( 79 ) 



( 80 ) 



( 81 ) 



where Tc. is the chamber temperature and ~tZ is the time measured 
from the Instant the particles enter the convergent portion of the noz- 
zle. The conductivities of the solid and gas are L s and re- 

pectively. 

For steady state, with spherical symmetry, the equation for 
heat conduction in the gas is 



b y~ 



( 



(~ 



i c>T 



— O 



( 82 ) 
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which has the solution 



T= A. +- B 

r 

At h = °o , T= f and thus 

T= which gives A = cl 



b=t 3 



Thus 



and at , 



fix') 

\ ^ r / r=a_ + o 





(83) 



For quasi-steady heat flow in particle and gas, substitute 
from equation into 81 to obtain 



^ T r T ^ <84) 

This will be used for the approximate boundary condition for the heat 
conduction equation of the particle. The solution to these equations 
(79 t 80 f 84 ) obtained by making use of the Laplace trans- 

form. 

Define 

f\e~ f F( t)Jr = FOp) (85) 

4> 



and 



V'-' 1 — <3? 



(86) 
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Then equation 79 becomes 

+ ^ To=o 

Since there is no source at f —<£> , \~z. = <=> at — ° and 

. f _/ °_ r _ _/ 

therefore ~7~ “• o where T ~ c r- — ~ 



(87) 



which means that 



t3 - 3 — > o 



as y~ 



( 88 ) 



The solution of equation 87 is 

r - f$r 

o — Ae + Be -hT c r 



(89) 



Substituting this into equation 88 



+ T„ 



* q r ' 

ho, ~t Be + Tc r 



(90) 



and to satisfy condition 83 , B — A t 



thus 



T = ®-' = Ixs'inUh 4 -T c 

V Y~ o 



(91) 



where ^ = \jfyP 

Substituting this into tlie boundary condition (equation 84 ), 
we obtain for the left hand side 

k s {z^r-^L r ]) 

J r -cv_ 
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and for the right hand side 



A 



Ol 



— ^5 F TV — S\ n k. a~ __ 

<3- 



-( T r T -)| ■ 



«- -T tt ] 



since r - t»- at >“ =• ■a— . Combining these equations, equa- 

tion 84 becomes 



ks- 



14 



■ fr CczLq a. Si n ^ f 



a- 



= V[t| - V 



(92) 



and solving for A 



A = ko ( Tu T c ) 



4 _ 

k* £ d£>s/v. ^ a. - S tt\k ^ -Jjj. j 

Substituting equation 93 into 91 gives 



(93) 



T = 



y ks * 1 



L ^ a. (t 

kj' 



+ I, 



(94) 



Let us now define an average sphere temperature T~ s as 



- CL. 



-r _ f mrSTit- 

»S — -r* ■ 



and hence 



-CL. 



71-^ r‘TJ, 



(95) 
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Then T s will be compared with the uniform sphere temperature in 
the heat balance equation for infinite (reference 1). 




Also define: 

B= 

- Si r\k ya. [\ - 

so that equation 94 may be written as 



T = B s±sX_$jt , Tc 
r 



(97) 



Then substituting and evaluating equation 95 




(98) 



Because large values of time T. correspond to small values 
of , expand the expression for 77 In a power series in <^_ to 
obtain the form correct for small , i. e. , 
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T r % 



Retaining only terms to order 




in equation 99 




which reduces to 

Tyr, - (77-7-4 1 (rf] (ioo) 

Recalling that ^ ^ ^ — ■ we find that for 

equal to infinity, equation 100 reduces to 

( 101 ) 

and talcing the inverse Laplace transformation of this leads to 






(102.) 



which is identical with the assumption used in the heat balance of a 
particle of infinite thermal conductivity, i. e. , equation 96 . 

Reference 2 gives values for the thermal conductivity of 
aluminum oxide and its variation with temperature. Measured values 
are given for temperatures up to 1600°C, and one more point (1 300°C) 
is presented from extrapolated data. The thermal conductivity de- 
creases with an increase in temperature up to 1400°C, and then be- 
gins to increase again. The value given for 1400°C (the minimum 

value of 4s for this range) is .0125 — or 

cm sec C 
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36. 25 



BTU 



m 
o. 



ft hr W F 

To find a reasonable value for the thermal conductivity of the 
gas, the data from the numerical example in reference 1 was used, 
that is, 

C = 0. 500 BTU/lb °F , 

P 

^ = 1. 5 x 10' 6 lb. sec/ft 2 , 

Pr = 0.74 . 

To be most conservative in evaluating ( a value of 

ju. was used that would give a high ^ • Since f-i varies con- 

siderably with temperature and the Prandtl number is almost constant 
over a large temperature range, the thermal conductivity varies with 
the temperature in the same manner as /U. . Over the temperature 
range of interest, a power law of the following form is satisfactory. 









afc. 



(103) 



The maximum L ^ would be obtained when = 1c. • thue 

k Sf 



= 1.412 BUL . Sjl. 

r J-i £ >_ 

rr U r 



Thus the value of ^y'kj in equation *^0 will always be 
greater than 25. 7 , and thus one may neglect 1/5 in comparison to 



1 



With this simplification in equation and substituting 






for 



% 
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V r ‘ 



= (h-x)l 



3 



k 



3 



f 




(104) 



Talcing the inverse of the Laplace transform yields 




(105) 



Without the second term on the right, this equation is identical with 
equation 101, and the last term is a correction term. To investigate 
the size of the last term, rewrite equation 102 a3 



d Jl = 

d r 




(106) 



and differentiate to obtain 



ilii. = ocl iB. - J r s 



At 



l 4- 

dx At 



(107) 



where oi — £ cu~ • Then equation 105 becomes 



L- (T 3 -T s )-Xp s Ca. 






(108) 



The condition that the approximation used in the heat balance 
equation of a particle, equation 96 , is valid is 



A fe« J * 



W JTt 

dx 



ils 

dr 



( 109 ) 
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For a first approximation, one can say that 






T and 



are of the same order, while still unevaluated. The coef- 

-4 

ficient of the numerator is 1. 621 x 10 based on the constants evalu- 
ated earlier. Since can be at most the same order 

as ‘‘Mr and will probably be lees, we can say that condition 109 
is satisfied. 

Since this condition is satisfied, equation 96 is a valid ap- 
proximation which would be in error by much less than one per cent. 
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LIST OF SYMBOLS 



cross section area 
specific heat 




Mach number 
Prandtl number 

ratio of absorption cross-sectional area/geometrical 
cross-sectional area 

gas constant 

Reynolds number 

temperature 

axial velocity component 
radial velocity component 



radius of solid particle 
radius of combustion chamber 
internal energy /unit mass 
correction factor 
enthalpy 

thermal conductivity 
Rosseland mean free path 
mass flow rate 

number of particles /unit volume 
pressure; variable of Laplace transform 
rate of heat transfer 
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% 

f 

t 

* 

oc 

t 

6 

X 

6 > 

A 4, 

? 

r 

r 

4 > 

r 




radial distance; burning rate 
time 

axial distance 
mass fraction of particle 
ratio of specific heats 
indicator of a small quantity 
pressure ratio 

non-dimensional radial distance 

angle ( radians ) 

viscosity 

non-dimensional axial distance 

density 

diffusivity 

Stefan- Boltzmann constant 
temperature ratio 

non-dimensional axial velocity component 
non-dimensional radial velocity component 



Subscripts 

j) burning surface 

C chamber 

d drag 

g gas phase 

convection 
£ with lag 
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f 

r 

s 

t 

V 

o 

1 



constant pressure 
radiation 
solid particle 
throat 

constant volume 
no lag 

with small lag 



